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ABSTRACT 

The discovery of OI atoms and CII ions in the upper atmosphere of HD 209458b, made with 
the Hubble Space Telescope Imaging Spectrograph (STIS) using the G140L grating, showed that 
these heavy species fill an area comparable to the planet’s Roche lobe. The derived ~10% transit 
absorption depths require super-thermal processes and/or supersolar abundances. From subsequent 
Cosmic Origins Spectrograph (COS) observations, CII absorption was reported with tentative velocity 
signatures, and absorption by Si III ions was also claimed in disagreement with a negative STIS G140L 
detection. Here, we revisit the COS dataset showing a severe limitation in the published results 
from having contrasted the in-transit spectrum against a stellar spectrum averaged from separate 
observations, at planetary phases 0.27, 0.72, and 0.49. We find variable stellar Si III and CII emissions 
that were significantly depressed not only during transit but also at phase 0.27 compared to phases 
0.72 and 0.49. Their respective off-transit 7.5 and 3.1% flux variations are large compared to their 
reported 8.2±1.4% and 7.8±1.3% transit absorptions. Significant variations also appear in the stellar 
line shapes, questioning reported velocity signatures. We furthermore present archive STIS G140M 
transit data consistent with no Si III absorption, with a negative result of 1.7±18.7 including ~15% 
variability. Silicon may still be present at lower ionization states, in parallel with the recent detection 
of extended magnesium, as Mg I atoms. In this frame, the firm detection of OI and C II implying solar 
or supersolar abundances contradicts the recent inference of potential x 20-125 subsolar metallicity 
for HD 209458b. 

Subject headings: planetary systems - stars: individual: HD 209458 - ultraviolet: stars - techniques: 
spectroscopic - planets and satellites: atmospheres 


1. INTRODUCTION 

UV detection of both hydrogen and heavy species on 
the extended atmospheres of hot exoplanets orbiting in¬ 
side 0.1 AU from their stars provides a useful tool for 
understanding the properties of these uppermost plan¬ 
etary layers, the energy input from the immense stel¬ 
lar X-ray and UV (XUV) radiation and the processes 
of interaction with the incident stellar wind magnetized 
plasma. Provided information on the velocity distribu¬ 
tion can be estimated in the transmission signature, the 
relative abundance of heavy species can be constrained 
to in turn reveal compositional properties that also per¬ 
tain to the lower atmosphere and nature of the exoplanet. 
UV studies of exoplanets can thus provide unique infor¬ 
mation for their characterization that is usually unac- 
cessible with IR and optical observations. Most of what 
we understand about hot-Jupiter upper atmospheres is 
based on far-ultraviolet (FUV) Hubble Space Telescope 
(HST) transit observations and extensive modeling of 
HD 209458b, which orbits a GOV star at 0.047 AU. How¬ 
ever, as we show in this paper, the full HST FUV obser¬ 
vational potential of this target has yet to be reached. 
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Using the Space Telescope Imaging Spectrograph (STIS) 
with the G140M medium-resolution (R ~ 10,000) grat¬ 
ing, a ~10% absorption in Lyman-a was discovered 
by HI atoms in an inflated atmo sphere reach ing the 
Roche lobe distance of ~ 2.9 R„ i (jVidal-Madiar et all 
I2003L[20081: rBen-Jaffell[2007ll2008ll . A Lyman-a absorp¬ 


tion with an enhanced blue-shifted component was first 
reported from the a verage of three med i um-re solution 
G140M tr ansits by IVidal-Madiar et all (1200311 . while 
[Ben-JafO ( 2008) . who included a fourth partial tran¬ 
sit, found a more symmetric average absorption in the 
red and blue wings. Existing theoretical models have 
tried to reproduce the HI observations assuming the 
solar XUV input where the atmospheric inflation re¬ 
sults from heating at the base of the thermosphere 
around 0.01-1 /ibar, yielding a mean temperature of order 
10,000 K that drives a hydrodynamic outflow. Most mod¬ 
els find enough thermal atomic hydrogen to reproduce 
the s ymme tric absorption profile originally r eported in 


Ben-J affel (120081) for the HD 209458b transits (iBen-Jaffell 
20081: Ben-Jaffel fe Hosseinil 120101 : iTrammell et all 120111 : 


Koskinen et all 120101 l2013al ). However, some issues 


remain unresolved regarding the exact total thermal 
H content that depends on the atmospheric solar or 
supersolar abundances assumed ([B en-Jaffel & Hosseinil 
1201(1: IKoskinen et aD [Mil I2013allbfl . on th e streng th 
of the outflow and the tidal effects fe.g.. lLammer et all 
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2001 lYeiai200i iTian et al.ll200l ICUircfa Munoz! |2007 


Murrav-Clav et alj|2009b Stone fc Progall2009UGuol 2011 


20131 : iBourrier et all 12014) ) . on the geometrical distribu¬ 
tion of the gas which can be global or quite likely confined 
at low (magnetic) lat itudes due to the presence o f a plan¬ 
etary magnetic field dT rammell et, al Jl20TTl 1201 4 1 Adams! 


etary magnetic neld (trammel! et al. z!Ui4; Adams 

2011b iKhodachenko et al.l l2012t iCohen fe Gloceil 120121 


Owen fc Adamsll2014lh as well as on the presence of non- 


thermal processes affect ing the energetics and line broad¬ 
ening (|Ben-.Taffel fc Hosseinil DOlTilhand on the s t rengt h 
of the stellar XUV heating " ([Ballester et all 1201 5). 
Char ge exchange with the i ncident stellar wind pro¬ 
tons (jHo lmstro m et al.l 120081 : iTremblin fc Chiand [20131 : 
IBourrier fc Lecavelier 2013T) a nd stellar radiation pres- 
sure (lYidal-Madiar et al.ll2003l ) are also key mechanisms 
that predict an asymmetric Lyman-a transit absorption 
profile, but the relevance of these proc esses depen ds on 
the inclusion of natural line br oade ning (lB en-Jaffell2008fl 
and self-shielding (IBourrier fc Lecavelied 120131 ), and a 
comprehensive analysis is require d including the presence 
of a planetary magnetic field dKislvakova et al.l 120141 ) 
and the inherent absorption by the interstellar medium 
thay may erase the asymmetric absorption feature if the 
Doppler shift is not strong enough. Probably the inter¬ 
stellar medium (ISM) absorption and the relatively low 
signal-to-noise ratio (S/N) limit the HI Lyman-a data 
analysis and interpretation, yet all studies agree on the 
important result that an area equivalent to the size of the 
Roche lobe is filled on HD 209458b, a finding that has 
many implications on our understanding of these atmo¬ 
spheres when con nected w ith transi t signatures in other 
species (iBen-Jaffel fc Hosseinil [20loi h However, we must 
stress here that available observations already show that 
we are close to making a final diagnostic on the very na¬ 
ture of the transit absorption as shown by the comparison 
between light curves derived for the blue and red wings 
of the HI Lyman-a line (e.g., Fig. 3 in iBen-Jaffelll2008h . 
For the rel atively good signal to noise reported thus far 
(|Ben-Jaffell l2008n . repeated new STIS G140M medium- 
resolution transit observations of HD 209458b should 
help converge our determination of the HI Lyman-a 
transit absorption profile for this planet, and definitely 
elucidate any asymmetric spectral absorption signature. 
Additionally, time variation may be present in the HI 
Lyman-a spectral signature of HD 209458b, as can be 
seen in Fig. 3b of (jBen-Jaffelll2007t b Similarly, time vari¬ 
ation has been clearly identified in the HI Ly man-a ab¬ 
sorpt i on of the hot-Jupiter HD 189733b dLecavelier et al.1 
12012 IBourrier et aHl20H . On HD 189733b, one STIS 
observation showed no HI Lyman-a absorption while on 
another date it showed absorption that was enhanced in 
the blue wing. A STIS observation of GJ 436b has also 
shown an extended hydrogen atmosphere on this warm 
Neptune, revealing a post-egress signature from a hy¬ 
drogen tail with juy enhanced blue-shifted HI Lyman-a 
absorption dKulow et al.ll2014f) . 

On HD 209458b, the STIS detection of extended neu¬ 
tral oxygen atoms and singly ionized carbon (the two 
dominant forms of the species for this environment) also 
extending to near Roche lobe distances has also been 
of great importance, since these components can pro¬ 
vide addi tion al critica l diag nosti cs of the up per at rno- 
sphere (lYidal-Madiar et ahl 2004 IBen-Jaffel fc Hosseinil 
l2010f) . Their detection shows that heavy species are ef¬ 


fectively entrained in the hydrodynamic outflow of the 
lighter hydrogen gas (through collisions with the neu¬ 
tral hydrogen ato ms and the protons) and overcome dif¬ 
fusive separ ation (jGarcfa Munod 120071 : iKoskinen et al.1 
I2010L [2fTl3alTbl) . The transit absorption depths in the OI 
1304 A and CII 1335 A multiplets derived from the STIS 
G140L observations of HD 20 9458b were, respectiv ely, 
12.8±4.5% and 7.5±3.5% (lYidal-Madiar et al.ll2004lb or 
10.5 ± 4.4% and 7.4 ±4.7% f rom a re-visit of the observa¬ 
tions dBen-.Iaffel fc Hosseinill20T0 l and correspond to 2.8 
and 2.2-sigma, or 2.4 and 1.6-sigma detections. These 
large absorption depths cannot be explained by 10,000 K 
thermal populations at solar abundances since line pro¬ 
files that would result from thermal and natural broad¬ 
ening for vertical distributions from 1-D hydrodynamic 
models are not wide enough to significantly absorb the 
broad stellar lines (with solar line widths) during transit. 
This is explored in various papers (see below). Since the 
line absorption profiles and thus the velocity distribu¬ 
tions of the absorbers were unresolved with STIS (using 
the G140L grating at R~1000), substantially different 
interpretations have been put forth. 

To broaden the absorption profile while keeping solar 
abundances, energetic super-thermal populations have 
been invoked (with effective temperatures x 5-100 larger 
than a ^10,000K background H-dominated gas), such as 
from processes that preferentially energize minor species 
involving chemistry, radiation, and wave-particle inter¬ 
actions, or from turbul ence a nd the stellar wind inter¬ 
action (IBen-Jaffel fc Hosseinil [20Toh . Ben-Jaffel & Hos- 
seini have shown that the final transit absorption pro¬ 
file is strongly dependent on the vertical distribution of 
the super-thermal population, a diagnostic that could 
be used if medium-resolution transit spectra are ob¬ 
tained. Alternatively, to enhance the optical thickness of 
the thermal population, some hydrodynamic models find 
that x3-5 supersolar abundances would be required ac¬ 
companied by a hotter thermosphere from tw ice t h e stel¬ 
lar XU V input than from the mean Sun (IKoskinen et al.1 
l2013bth (A previous hydrostatic approximation invoked 
quite large x4-40 supersolar abundances with an unrea¬ 
sonably stron g XUV flu x x5-10 higher than the mean 
Sun: IKoskinen et al.ll20ldh . However, recent findings are 
now showing an unexpected long-term low activity on 
the star HD 209458, and this low activity does not sup¬ 
port scenarios that invoke even moderately high XUV 
fluxes, such as those of solar mean or maximum activity, 
to heat and inflate t he atmosphere to the e xtent required 
by the observations dBallester et al J 1201511 . As there ap¬ 
pears to be an energy crisis in the upper atmosphere 
of this planet due to the low XUV input, new informa¬ 
tion on super-thermal processes that may operate in the 
upper atmosphere of this planet are needed, and such 
evidence is uniquely accessible with transit observations 
that resolve the line absorption profiles and thus velocity 
and energy distribution of the heavy species. 

The Cosmic Origins Spectrograph (COS) on HST pro¬ 
vides higher sensitivity than STIS, and it can improve the 
S/N of the transit absorption in CII 1335A and in OI 
1304 A as well as resolve the multiplets and the line ab¬ 
sorption profiles at medium spectral resolution (with the 
G130M grating at R~20,000) as has been demonst rated 
for the case of HD 189733b (IBen-Jaffel fc Ballester!l2013l . 
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l2014li . FUV observations of the HD 209458 sy stem were 
made with COS G130M bv lLinskv et al.1 (1 20 1(1 ). Absorp¬ 
tion by the planet during transit was reported in the CII 
1335 A doublet lines with a depth of 7.8 ± 1.3% in the 
-50 to +50kms _1 interval, similar to the line-integrated 
absorption detected with STIS. Tentative velocity signa¬ 
tures were reported, with absorption at -10 to lSkms^ 1 
and less certain features at -40 and +30-70 km s -1 , 
and no absorption at zero velocity. The larger veloci¬ 
ties would agree with the presence of an upper atmo¬ 
spheric hot populati on, or of plasma wi t h large motions 
around the planet dLinskv et all 120101 : [Koskinc n et al.l 
I2013bf) . with great implications about the properties of 
HD 209458b’s upper atmosphere;_ 

From the same COS dataset, iLinskv et al.l (120100 also 
reported an 8.2 ± 1.4% transit absorption by Si III ions 
in the 1206.5 A line. In sharp contrast, no Si III 1206.5 A 
absorption was detected in the average of the four HST 
transit observations made with STIS G140L that de¬ 
tected OI and CII, where a negative result of O.OIj/q% 
was derived in Si III dVidal-Madiar et aLll2004D . To fit 
the COS Si III results. IKoskinen et al.l (|2013bi ) found sim¬ 
ilar x3-5 supersolar abundances and higher stellar en¬ 
ergy input requirements for the upper atmosphere, sim¬ 
ilar to the inferences made for the OI and CII. 

The composition at high atmospheric altitudes (as 
modified by the stellar ionization, thermal balance, and 
charge-exchange processes) reflects the species and abun¬ 
dances at the base of the thermosphere, and it is 
thus diagnostic of the properties of the lower atmo¬ 
sphere an d overall planetary abundances and energy bu d¬ 
get (e.g.. iGarcfa Munod [20071 : IKoskinen et al.1 [2013al lbl: 
lLavvas et al.l 120141) . This is particularly the case for a 
hot Jupiter with its strong 3D atmospheric circulation 
and l arge effective eddy mi xing (e.g., IShowman et al.l 
120091: iParmentier et al.H2013ll . and where the detection 
of extended OI and CII on HD 209458b itself shows 
that diffusive separation of the heavy species in the up¬ 
per atmosphere is impe ded by the hydrodynamic outflow 
(IKoskinen et al.|[20 1 3 al ibi) . For OI and CII the key lower 
atmospheric components are H 2 0, that does not con¬ 
dense on hot Jupiters, and CO, the domin ant carbon - 
bearing species on hot Jupiters (e.g.. IMoses et al.ll2011l) . 
For refractory species like silicon, the situation is more 
complex and unce rtain. For silicon to b e present in 
the thermosphere, IKoskinen et al.l (l2013bl) find that it 
cannot condense into silicate clouds in the lower atmo¬ 
sphere, and that the vertical transport should be strong 
to keep any condensate particles lofted. Furthermore, sil¬ 
icon must be predominantly in the Si II ionization state, 
and charge exchange with protons (and He + ) can then 
supply a significant Si III population to be detected. It 
may well be that silicon is not effectively removed by 
vertical and horizontal cold-traps on this planet, given 
t he recent d etection of JVlgl in the upper atmosphere 
dVidal-Madiar et ahll2013D , but we must first address the 
validity of the reported COS detection of silicon ions, at 
least as Si III, given the contradiction with the negative 
STIS G140L results. 

In this paper we re-visit the COS observations of the 
HD 209458 system given that these observations were not 
made in the standard fashion of observing immediately 
before and/or after the transit to contrast the in-transit 


fluxes with the actual stellar fluxes around the time of 
transit. Instead, the in-transit stellar fluxes were ratioed 
against the average fluxes measured on three separate 
HST visits that were obtained over a period of one to two 
weeks with the planet off transit. The limitation of this 
approach is that the stellar activity, even if it is low, can 
cause significant flux variations in the FUV line emissions 
such as over a stellar r otati on peri o d, and fro m sto chastic 
or flare activity (e.g.. iParke Lovd fe Franck I2014T ). This 
is indeed the case for the Sun for which a variation of 
~ 27% was reported for the Si III 120 6.5 A line during 
a rotation period (e.g., Table 3 in lBen-Jaffel fc Ballesterl 
I2013D . It is thus possible that the transit absorptions 
derived from the COS data could be related to stellar 
variation rather t ha n to the true planetary transit signa¬ 
ture. ILinskv et al.l ( 2010 ) deduced their approach to be 
valid because the flux level and line shape in the Si IV 
1393.7 A line was similar in the in-transit data and in the 
averaged off-transit data. We find, however, that this use 
of the averaged off-transit data requires reassessment. 

In Sec. 2 we re-evaluate the COS observations in the 
Si III 1206.5 A line, CII 1335 A doublet, Si IV 1393.7 A 
line, and 1364.8-1390.7 and 1405.6-1423.6A continua, 
taking advantage of the temporal information of the 
data, both from the time-tagging and from comparison 
of the average stellar fluxes per HST visit. We find in¬ 
deed significant time variation in the stellar emissions 
that render the previously derived transit absorptions 
for Si III 1206.5 A line and the CII 1335 A doublet from 
this COS dataset unreliable. Complications in the usage 
of the COS dataset for a derivation of velocity signa¬ 
tures in the absorption are also addressed. Moreover, we 
also present previously unpublished Si III transit obser¬ 
vations of HD 209458b from a re-visit of the original STIS 
G140M dataset that discovered extended HI in Lyman-a 
absorption, since these exposures also sampled the Si III 
1206.5 A line. We find no evidence of a Si III transit ab¬ 
sorption, in agreement with the first STIS G140L find¬ 
ings. We present and discuss the STIS data in Sec. 3. 
We also discuss some issues related to the conditions that 
would be needed in the lower atmosphere to obtain silicon 
ions in the upper atmosphere. We finish by establishing 
the need for COS FUV transit observations of this ex¬ 
oplanet, in particular for the unexploited diagnostics in 
the OI 1304 A and CII 1335 A transit signatures. 

2. COS OBSERVATIONS 
2.1. Observations and data reduction 

HD 209458 was observed with COS on Sep-Oct 2009 
during four HST visits, two or three HST orbits each, 
with the planet around orbital phases 0=0.27, 0.72, 0.0, 
and 0.49 (Table 1). The COS G130M grating was used 
covering the ~1130-1460A region at medium spectral 
resolution (R~20,000 ~15km/s, with 0.01 A sampling). 
Throughout each visit the exposures were made alter¬ 
nating between four grating central wavelength settings 
referred to as 1291, 1300, 1309, and 1318. COS FUV 
spectra are sampled by two detector segments, A and B, 
that are separated by a small gap. Some of the grating 
settings either did not sample the important 011304 A 
triplet or had the lines near the edge of the detector (Ta¬ 
ble 1), so this stellar emission could not be evaluated 
as was possible for transit observations of HD 189733b 
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TABLE 1 

Summary of COS observations 


Exp. 

No. 

Start time 
in 2009 
(UT) 

Exp. 

time 

(sec) 

Planet 

mid. 

orbital 

phase 

Grating 

wavelength 

setting® 

01 

Sep 

19 

10:10 

2340 

0.255 

1291 

02 

Sep 

19 

11:36 

0955 

0.270 

1300 

03 

Sep 

19 

11:55 

1851 

0.275 

1309 

04 

Sep 

19 

13:11 

0560 

0.288 

1300 

05 

Sep 

19 

13:24 

2235 

0.294 

1318 

06 

Sep 

24 

13:12 

2340 

0.710 

1291 

07 

Sep 

24 

14:38 

0945 

0.724 

1300 

08 

Sep 

24 

14:57 

1787 

0.730 

1309 

09 

Sep 

24 

16:14 

0925 

0.743 

1300 

10 

Sep 

24 

16:33 

1798 

0.748 

1319 

11 

Oct 

02 

14:32 

1045 

0.993 

1291 

12 

Oct 

02 

14:53 

1096 

0.997 

1300 

13 

Oct 

02 

15:55 

1400 

0.010 

1309 

14 

Oct 

02 

16:22 

1405 

0.015 

1319 

15 

Oct 

18 

10:51 

1057 

0.490 

1291 

16 

Oct 

18 

11:18 

1093 

0.494 

1300 

17 

Oct 

18 

12:20 

1401 

0.507 

1309 

18 

Oct 

18 

12:46 

1405 

0.512 

1319 


a The spectral ranges covered per grating wavelength set¬ 
ting respectively for detector Segments B &; A are: (1291) 
1132-1274 & 1291-1433 A; (1300) 1141-1283 &1300-1442A; 
(1309) 1153-1294 & 1309-1449 A; and (1318) 1163-1303 
&1319-1459 A. 

(]Ben-Jaffel fe Ballesterlf2013L I2014D . Note also that the 
CII 1335 A feature is a triplet, but it is detected as a 
doublet with a line at 1334.5 A and two unresolved lines 
at 1335.7 A. 

The time-tagged exposures were reprocessed into 100 s 
segments, but since the S/N was low, the segments were 
averaged into 1000 s sub-exposures so that the temporal 
flux variations (depicted in Fig. 1) were not dominated 
by the photon noise level. For the extraction of the stellar 
emission line fluxes, a weak stellar continuum was first 
subtracted from each spectrum based on a linear fit to 
the spectral regions devoid of obvious emission l ines iden- 
tified in a high-resolution solar spectrum (iCurdt et al.1 
l200l . For the Si III 1206.5 A line we also subtracted 
a linear fit of the blue wing of the stellar HI Lyman- 
a line which is very broad as detected by COS. This 
is a relatively small correction compared to the poten¬ 
tial planetary absorption, and the associated errors have 
been included in the propagated errors. The stellar line 
fluxes, depicted in Fig. 1, were extracted by co-aligning 
the given line in all the spectra, and integrating the spec¬ 
tral f luxes out to abo ut 10% of the averaged peak inten¬ 
sity (B en-JaffeII120071) . This integration spanned about 
+/-28km/s. Photon noise errors were propagated in all 
reduction steps and problematic pixels (with dq_wgt=0) 
were ignored. 

2.2. Absolute Wavelength Scale 

The determination of a common absolute wavelength 
scale, and thus of any planetary velocity signature re¬ 
ported from the COS dataset, was complicated because 
of having alternated between the four G130M central 


wavelength settings during each visit. The original 
purpose of that setting was to minimiz e fixed-patte rn 
noise and grid-wir e shadow effects (Fr ance et al.ll2010t 
iLinskv et al.1 [20 10b while for most other transit obser¬ 
vations of exoplanets keeping instrumental settings fixed 
has proven best for measuring relative changes due to 
the planetary transmission alone. The grid-wire shadow 
effects are now automatically corrected by the COS 
pipeline software, as in the CALCOS 2.14.4 and 2.18.5 
pipeline versions used for the time-tagged flux analysis 
and the wavelength study, respectively. The main com¬ 
plication is that the repositioning of the gratin g is not 
exact ly repeatable due to thermal flexures dShaw et al.1 
120091) . Although in principle the absolute wavelength 
scale for each exposure could be set by co-aligning the 
CII 1334.5A line based on the dip from the ISM ab¬ 
sorption of this line as done bv ILinskv et al.1 (12010I) . we 
found that the S/N is too low to accurately find the cen¬ 
ter of this absorption in each exposure. To co-align the 
spectra, for the Segment-A data we compared co-aligning 
the spectra against each line of the CII 1335 A doublet 
and also against the full doublet. Significant differences 
were found when using only one line or the full doublet, 
when the data were smoothed or not, and when differ¬ 
ent exposures were chosen as the template. Comparing 
results from fitting against exposures 1 and 15 (Table 
1), using only the CII 1334.5 A line or the full doublet, 
and with or without smoothing, we found offsets in the 
centering of the ISM 1334.5 A line absorption of about 
lOkms -1 (about 5 pixels). This error, though somewhat 
smaller than the 15kms _1 spectral resolution, produces 
large differences in the resulting line profiles and in the 
corresponding transit absorption line profile. In Fig. 2 
(Sec. 2.4) we show a sample result of line profiles from 
co-alignment against exposure 15, using the full doublet, 
and with the data box-car smoothed by 5. 

The Si III 1206.5 A line data, sampled by the detector 
Segment B, was co-aligned separately. The wavelength 
shift found from the CII 1335 A doublet in the Segment 
A data differed from that found for Segment B by up to 
±4 pixels (while Linsky et al. 2010 applied the same shift 
to data from both segments). This turned out to be ex¬ 
pected since the COS pipeline calibration calculates the 
spectral offset for each detector segment independently. 
The SHIFT1A and SHIFT1B in the processed data differ 
by about 4-7 pixels, similar to the offsets that we found. 

2.3. Time variation of the stellar fluxes 

Figure 1 shows the time-resolved fluxes for the Si III 
1206.5A line, CII 1335A doublet, SiIV 1393.7A line, 
and 1364.8-1390.7 and 1405.6-1423.6 A continua. Black 
symbols show the fluxes from the 1000 s sub-exposures, 
and gray symbols show the averaged fluxes per HST visit. 
For each feature, the fluxes have been normalized by 
the corresponding average value for visit 4 (</> ~0.49) 
to ease the interpretation of the results. The error bars 
for the sub-exposures are propagated statistical photon 
noise, while the error bars for the visit-averaged fluxes 
include the non-statistical flux scattering within a given 
visit defined as the standard deviation from the mean. 
The variations within an HST visit can be stel l ar but 
there can also be instrumental effects (e.g.. ILinskv et ahl 
12013) . Since the intra-orbit flux variations do not seem 
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0 1 x 1 0 4 2x 10 4 3x10 4 

Observation Time (s) 


Fig. 1.— Normalized fluxes for the 1000 s sub-exposures versus 
observation time for the stellar emission features Si III 1206.5 A 
line, C II 1335 A doublet, SiIV 1393.7 A line, and the 1364.8-1390.7 
& 1405.6-1423.6 A continuum. The four HST visits are labeled on 
the top panel by the planetary orbital phase. For each visit the 
times have been offset from their arbitrary zero-time by 0, lxlO 4 , 
2xl0 4 , and 3xl0 4 s, respectively. For each emission feature, the 
average flux per visit is over-plotted as a gray circle at the middle 
of the visit. 

repeatable from orbit to orbit, they cannot be systemat¬ 
ically corrected for. The stellar fluxes averaged per visit 
are provided in Table 2 with both the photon-noise errors 
as well as errors that include the intra-visit scattering. In 
the discussion below we refer to the photon- noise errors 
for direct comparison with the results by iLinskv et al.1 
« ■ 

To evaluate the validity of deriving a transit absorp¬ 
tion from this COS dataset we look at the normalized 
fluxes averaged per HST visit (Table 2). For the CII 
1335 A doublet we find that although the flux during 
transit {4> = 0.0) is low at 0.880±0.021, the flux is also 
relatively low for cj) = 0.27 at 0.929±0.017. In contrast, 
the flux for (f) = 0.72 of 0.982±0.018 is similar to the ref¬ 
erence value of 1.0±0.023 for (j) = 0.49 within the errors. 
Although the transit depth relative to the mean of the 
off-transit fluxes (at cj> = 0.27, 0.72, and 0.49) would be 
9.3±2.4%, which is some what deeper than the 7.8±1.3% 
reported bv lLinskv et al.1 (l2010f) . the range of depths rela- 


TABLE 2 

Normalized stellar fluxes from COS dataset, averaged 
per HST visit 


Spectral 

feature 

Planet 

mid. 

orb. 

phase 

Norm. 

flux® 

Photon 

noise 

error 

Intra¬ 

visit 

std 

error b 

Photon 

and 

intra¬ 

visit 

error b 

Si III 1206.5 A 

<t> = 

0.27 

0.837 

0.018 

0.054 

0.057 

line 

<t> = 

0.72 

0.964 

0.020 

0.056 

0.059 


<t> = 

0.00 

0.853 

0.023 

0.079 

0.082 


<t> = 

0.49 

1.000 

0.025 

0.074 

0.078 

CII 1335 A 

<t> = 

0.27 

0.929 

0.017 

0.060 

0.063 

doublet 

<t> = 

0.72 

0.982 

0.018 

0.108 

0.109 


<t> = 

0.00 

0.880 

0.021 

0.065 

0.069 


cj, = 

0.49 

1.000 

0.023 

0.064 

0.068 

Si IV 1393.7 A 

<t> = 

0.27 

0.943 

0.029 

0.125 

0.128 

line 

0 = 

0.72 

0.933 

0.029 

0.070 

0.075 


b = 

0.00 

0.929 

0.036 

0.089 

0.096 


b = 

0.49 

1.000 

0.038 

0.131 

0.136 

1364.8-1390.7 & 

cf> = 

0.27 

1.032 

0.015 

0.063 

0.064 

1405.6-1432.6 A 

</> = 

0.72 

1.008 

0.015 

0.069 

0.070 

continuum 

<f> = 

0.00 

0.965 

0.019 

0.058 

0.061 


</> = 

0.49 

1.000 

0.019 

0.083 

0.085 


“ The fluxes have been normalized with respect to the av¬ 
erage of visit 4 (<f> = 0.49), given by: 2.12 X 10 -15 for 
Si III; 2.33 X lO" 15 for CII; 0.96 X 10“ 15 for Si IV, and; 
4.87 X 10 —15 ergs cm -2 s —1 for the continuum. 

b The intra-visit error is the standard deviation of the scat¬ 
tering of the stellar flux during the given visit. 

five to the different average fluxes for the off-transit visits 
spans from 5.3±2.9% to 12.0±2.9%. The detection and 
depth of the CII 1335 A transit signature of HD 209458b 
associated thus far with the COS instrument remains un¬ 
confirmed given the inappropriate observing method that 
was employed. 

The Si III 1206.5 A line also shows a strong time vari¬ 
ation. A low flux of 0.837±0.018 is observed at cj) = 0.27 
and a similar value of 0.853±0.023 is seen during transit. 
Furthermore, the flux of 0.964±0.020 observed at <f> = 
0.72 is comparable to the reference value of 1.00±0.025 
at (j) = 0.49 within the errors. Deriving a transit ab¬ 
sorption for Si III is therefore misleading, given the large 
time variations. A detection of Si ++ ions in the upper 
atmosphere of HD 209458b is not valid from this COS 
dataset. 

The 1364.8-1390.7 & 1405.6-1432.6 A continuum con¬ 
sists of the photospheric continuum (that decreases 
towards shorter wavelengths) and a white continuum 
emission from the chromosph ere (that increases to¬ 
wards shorter wavelengths: ILinskv et ahll2012D . No large 
continuum-flux variations are observed in tandem with 
the line emissions, indicating that no large common¬ 
mode instrumental effect was at play. This chromo¬ 
spheric component is expected to vary with stellar ac¬ 
tivity, yet the variability at these wavelengths shoul d be 
smaller than in any of the emission lines ('e.g. JSnow et all 
f201fT ). Although in the observation at (j> = 0.27 the con¬ 
tinuum did not appear relatively low as in the other fea¬ 
tures but was instead somewhat larger, and during the 
transit visit it was relatively low as could be due to the 
transit of the planetary disk, the fluxes are all similar 
within the errors. 

For the Si IV 1393.7 A line (significantly brighter than 
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TABLE 3 

HD 209458 average stellar fluxes from COS dataset, and 

APPARENT VISIT-TO-VISIT VARIATIONS FROM OFF-TRANSIT DATA 


Spectral 

feature 

Flux' 1 ’ 6 

(io - 15 

ergs/cm 2 /s) 

%Variation 6 
as a 

std. dev. 

%Variation 6 
as a 

max/min 

Si III 1206.5 A line 

1.98±0.15 

7.5% 

19.4±8.3% 

C II 1335 A doublet 

2.26=1=0.07 

3.1% 

7.7±5.4% 

Si IV 1393.7 A line 

1.471=0.05 

3.1% 

7.2=1=10.7% 

1364.8-1390.7 & 
1405.6-1432.6 A cont. 

4.94=1=0.07 

1.4% 

3.2=1=5.0% 


" Mean spectral flux and standard deviation. 


° Errors include intra-visit noise. 

the other 1402.8 A line of the doublet), we find that the 
fluxes are relatively low not only at both <j> =0.0 and 
0.27, but also at <j> = 0.72 unlike for Si III and CII for 
which the fluxes at 0 = 0.72 were higher and comparable 
to the fluxes at <t> = 0.49. Thus, the stellar Si IV emission 
turns out not to be an accurate proxy for the CII and 
Si III emissions. 

The variations of the emissions as observed in the off- 
transit visits are listed in Table 3. Clearly, the 7.5% 
variation (or 19.4±8.3% max/min) of the Si III 1206.5 A 
flux is quite large, more so than for the CII 1335A or 
Si IV 1393.7 A fluxes although there may be some overlap 
if the intra-visit scatter is included in the errors. 

It is interesting to note that the 3.1% variation (or 
7.2±10.7% max/min) of the Si IV 1393.7 A line is smaller 
than the 7.5% variation of the Si III 1206.5 A line 
(or 19.4±8.3% max/min) for the same set of observa¬ 
tions (although again they overlap within the larger 
error bars). Since the SiIV 1393.7A emission origi¬ 
nates from a hotter layer in the stellar atmosphere (at 
logT = 4.75 for the Sun) than the CII 1335 A and Si III 
1206.5 A emissions (at logT = 4.10 and 4.25, respec¬ 
tively; |WoodIIeFaI][2()0(J it would be expected in gen¬ 
eral to vary more strongly or at least comparably to the 
Si III 1206.5 A line (based on more variability at higher 
temperature regions). 

Solar observations indicate that the Si IV 1398 A dou¬ 
blet variability is comparable in general to that of the 
Si III 1206.5 A line, or somewhat larger at high ac¬ 
tivity as known from short-term and lon g-term obser¬ 
vations and from observ ations of flares (Woo ds et al.l 
120001: iBrekke et al.l 119961b However, a recent evalua¬ 
tion of the solar variability in the FUV emission lines 
rel evant t o ho t-Ju p iter t r ansit studies has been made 
by iBen-Jaffel fe Ballester! (|2013T ) using 2003-2007 spec¬ 
tra from the Solar Stellar Irradiance Experiment (SOL¬ 
STICE) instrument on-board the Solar Radiation and 
Climate Experiment (SORCE) obtaining more detailed 
or relevant results for exoplanet FUV HST work. The 
data were divided into periods of low, medium, high, and 
extreme flare activity, where the lower activity should 
include the emission from plages and enhanced network. 
In three out of four flare level categories, the variabil¬ 
ity of the Si IV 1398 A doublet is comparable to that 
of the Si III 1206.5 A line, and somewhat larger for the 
cases o f high activity in agreement with the findings by 
IWoods et all (l2000ih In that study, there was a sepa- 


rate evaluation for t he S i IV 1393.7 A line (see bottom 
part of Table 3 in IBen-Jaffel fc Ballester! I2013D . The 
solar Si IV 1393.7A line showed, in contrast, that the 
variability with solar 27-day rotation and 11-year cycle 
was somewhat smaller, at 22 and 60%, respectively, com¬ 
pared to 27 and 73% for the Si III 1206.5 A line. (These 
variabilities are also standard deviations from the mean.) 
The solar lines do show more variability that the COS 
results for HD 209458, but the COS dataset is extremely 
limited to three samples. Here we present results for the 
Si IV 1393.7 A line and not the do ublet, since this was the 
line used by iLinskv et ali (|2010f ) as a proxy for the stel¬ 
lar activity given that it is about twice as bright as the 
second line in the doublet. In a separate analysis of the 
COS dataset (not presented here) in which we have not 
divided the exposures into time-tagged segments, we see 
that the variation in the full doublet is larger and more 
comparable to that of the Si III 1206.5 A line, due to in¬ 
clusion of the m ore variable and/or nois ier second line 
in the doublet. iParke Lovd fc France! (120141) report on 
temporal variations on FUV fluxes of sunlike stars in 60- 
second time intervals. For HD 209458 they report upper 
limits for the variation in the Si III 1206.5 A line and Si IV 
1398A doublet emissions (which is not surprising given 
the relatively low S/N of the data), and their upper limit 
variation for the Si IV doublet is a bit larger than for the 
Si III line. Therefore, from the various lines of evidence 
presented, it is unclear at this point if the relatively low 
visit-to-visit variation of the Si IV 1393.7 A line compared 
to the larger variation in the Si III 1206.5 A line is signifi¬ 
cant and distinct for this star compared to the Sun. The 
degree of intra-visit scatter for both features seems com¬ 
parable, however, although the COS data is extremely 
limited. Future observations should shed more light into 
the simultaneous temporal FUV flux variations on this 
star. 

Finally, it is also interesting that the Si III 1206.5 A and 
CII 1335 A emissions may have shown a stellar rotational 
variation. Visits 1, 2, 3, and 4 were made on days 0 
(reference date), 5, 13, and 29. Using the stellar rotation 
periods of ~ 11.2 - 14.1 days derived from the Rossiter- 
McLaughlin effect an d the stellar line wid ths provided 
in the literature (c.f.. IB allester et al.ll2015lh these visits 
corresponded to stellar rotational phases of 0 (reference 
phase), 131 T65, 336-63, and 20-212°. Therefore, visits 
1 and 3 that showed the lowest fluxes should correspond 
to similar stellar phases, around the reference phase 0°. 

We emphasize, however, that the COS dataset at hand 
is extremely limited, consisting of only a single in-transit 
observation visit and three off-transit visits on separate 
dates. To definitely ascribe a stellar rotational modula¬ 
tion to the Si III 1206.5 A and CII 1335 A stellar fluxes 
would be premature, since there can be stochastic ef¬ 
fects. Furthermore, to derive a reliable transit absorp¬ 
tion based on the observing technique used for the COS 
dataset, one would have to observe many times with 
the planet both in-transit and off-transit, to hopefully 
average-out potential rotational, stochastic, small flaring 
activity, and long-timescale magnetic cycle variations in 
the stellar fluxes. Such extensive observations are not 
viable with HST. Instead, the standard transit method 
of directly measuring the actual stellar flux immediately 
before and/or after transit is by far the most direct and 
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Fig. 2.— Sample shapes of the CII 1335.7 A line derived from the COS data of the four HST visits, box-car smoothed by 5, with error 
bars. The gray profile is the average off-transit data (0 = 0.27, 0.49, 0.72). To correct for non-repeatable wavelength offsets from alternating 
the grating setting, the data were co-aligned against exposure #15, using both lines of the CII 1335 A doublet. Visit-to-visit variation in 
the stellar line shape is apparent, but observations with a single grating setting should be more suitable to determine the stellar behavior. 


least ambiguous FUV transit observational method, as 
well as the one that would utilize the least HST resources. 
Already, the time variation found in the stellar fluxes in 
this COS dataset clearly demonstrates this argument. 

2.4. Problems with the reported velocity distribution of 
the absorbers 

COS G130M transit observations of HD 209458b have 
high enough sensitivity and spectral resolution to poten¬ 
tially determine the velocity distribution of the absorp¬ 
tion by the heavy species on the planet at similar S/N as 
in the STIS G1 40M observations of the HI Lym an-a line 
absorption (iVidal-Madiar et all I2003L 120081 : [B en-Jaffell 
120071 l2008t iBen-Jaffel fe Hosseinill2010l) . This has been 
demonstrated for both the CII 1335 A doublet as well as 
t he O 130 4 A tr ipl et transit o bservations of HD 189733b 
dBen-.Taffel fc Ballester! 1201 3b Figure 2 shows the av¬ 
erage CII 1335.7A line profiles for the four HST /COS 
visits of HD 209458. Although there may be real plane¬ 
tary atmospheric absorption during transit (<j> = 0), we 
find an apparent visit-to-visit variation in the average off- 
transit stellar line shapes that precludes a proper evalua¬ 
tion of the velocity signature in the planetary absorption. 
For these reasons of potential intrinsic time variation in 
the stellar line shapes, t he velocity dist ribution of the ab¬ 
sorbers reported by iLinskv et al.l (|2010tl for HD 209458b 
from the COS dataset at hand is invalid. We note that 
the apparent stellar line-shape variation may be related 
instead to errors in the co-alignment of the spectra that 
was needed since the multiple grating wavelength settings 
were not exactly repeatable as described in Sec. 2.2, but 
this cannot be fully explored with this dataset. 


TABLE 4 

Summary of STIS GI40M observations 


Exp. 

No. 

Start time 
in 2001 
(UT) 

Exp. 

time 

(sec) 

Planet 

mid. 

orbital phase 

o6e201010 

Sep 07 20:35 

1780 

0.980 

o6e201020 

Sep 07 22:13 

2100 

0.000 

o6e201030 

Sep 07 23:49 

2100 

0.019 

o6e202010 

Sep 14 21:12 

1780 

0.973 

o6e202020 

Sep 14 22:42 

2100 

0.992 

o6e202030 

Sep 15 00:18 

2100 

0.011 

o6e203010 

Oct 20 02:50 

1780 

0.969 

o6e203020 

Oct 20 04:20 

2100 

0.988 

o6e203030 

Oct 20 05:56 

2100 

0.007 


Proper COS transit observations are needed, with con¬ 
secutive data with the planet in and out of transit, and 
using a single grating wavelength setting (Sec. 2.2). 

3. STIS G140M AND G140L OBSERVATIONS: NO SI III 
DETECTION 

The HST STIS medium-resolution G140M transit 
observations that discovered t he extended hyd roge n 
atmosphere on HD 209458b (IVidal-Madiar et all 120031 
IBen-Jaffel 120071) also sampled the Si III 1206.5 A line. 
With a resolving power of R ~ 10,000, the line is spec¬ 
trally resolved and well separated from the wings of the 
Lyman-a line although the S/N is relatively low. Three 
transit observations were made (Table 4). Figure 3 shows 
the line fluxes per HST visit, integrated at 1206.5.07- 
1207.14 A and normalized to the average of the first two 
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Fig. 3.— Stellar fluxes in the Si III 1206.5 A line for three transits of HD 209458b observed with STIS G140M, the same dataset that 
discovered an extended HI atmosphere on this planet (IVidal-Madiar et al.H2003l ; IBen-Jaffell[20071 ). Diamonds for 7-8 Sep. 2001; triangles 
for 14-15 Sep. 2001; and circles for 20 Oct. 2001. The fluxes are normalized to the 1.56x 10 —15 ergs cm -2 s —1 average of the first two 
points. Th e data are integr ated per HST orbit (3 orbits per transit). The transit light curve is drawn for an atmosphere obscuring 8.9% 
of the star <IBen-Jaffelll200^T) . 


points of 1.56xl0 _15 ergscm~ 2 s -1 . No Si III transit ab¬ 
sorption is detected. The ratio of the average of the three 
points within 3000 s from mid-transit to the first three 
off-transit points yields a negative result of 1.7±18.7% 
absorption. The 1.7% absorption, if real, would repre¬ 
sent the 1.5% obscuration by the disk of the planet. 

Although the S/N is low, large temporal variations are 
clearly present, from visit to visit, and intra-visit. The 
std. dev. is 13% for the off-transit points, and 15% for all 
the data. Thus, as with the COS observations discussed 
in Sec. 2.3, we find significant variation in the stellar 
Si III emission. The 13% std. dev. variation of the off- 
transit data found with STIS is about twice the 7.5% 
value found with COS. 

The original STIS G140L observations sampled four 
partial transits and a negative detection in Si III 1206.5 A 
was reported with a transit absorption of 0.0t o 'o% by 
IVidal-Madiar et all (120041) . The non- detec tion was later 
confirmed bv iBen-Jaffel &; Hosseinil (l2010l) . Although a 
limitation of the low-resolution G140L data is the con¬ 
tamination of the Si III line by the blue wing of the stel¬ 
lar Lyman-a line, and the S/N in the (lower sensitivity) 
G140M data is low, taken together, there are now seven 
transits observed with STIS that do not show Si III ab¬ 
sorption. 

A summary of the current transit results of heavy 
species on HD 209458b with STIS and COS is provided 
in Table 5, including the ne w resu lts from STIS. We also 
include transit depths from lLinskv et all (l2010ll from the 
2009 COS dataset and those discussed in Sec. 2 for com¬ 
parison, although all of these depths are not valid since 


the stellar flux around the time of transit is unknown. 
For valid COS results new observations are required sam¬ 
pling the transit light curve. 

4. DISCUSSION: NEW LIGHT ON CONTRIBUTION OF UV 
STUDIES TO EXOPLANETARY ATMOSPHERIC 
SCIENCES 

The abundance of heavy species in the upper atmo¬ 
sphere has major implications on all pressure levels from 
the top to the bottom of a planetary atmosphere. The 
link of the FUV transit observations with properties of 
the planet and its lower atmosphere are often under¬ 
played. The mere detection of extended OI and CII 
absorption in HD 209458b’s upper atmosphere with STIS 
G140L (at low spectral resolution) already requires solar 
or supersolar abundances, a firm result that is in sharp 
disagreement with the recent infere nce of poten tial x 20 - 
125 low metallicity on this planet dMadhusudhan et al.1 
120141 ). This low metallicity has been raised as one of 
two explanations for an inferred low H 2 O abundance in 
the planet’s lower atmosphere: either a low metallicity 
and C/O < 1 since water would be the major oxygen¬ 
bearing species in this case, or; a high metallicity but 
C/O > 1. A X20-125 low metallicity would be highly 
unexpected on a hot Jupiter. (Note that HD 209458 has 
a solar metallicity [Fe/H1=- 0.0±0. 2; Mazeh et al. 2000.) 
The goal of iMadhusudhan et al.l 1 201 4!) is to explain a 
low contrast 1.4 /jm band absorption reported for this 
planet from HST Wide-Field Camera 3 transit data, yet, 
it is possible that s ome level of haze scattering is present 
(iDeming et al.ll20T3h . as well as that the data need fur¬ 
ther work (T. Evans et ah, in progress). Since species in 
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TABLE 5 

Current transit results for heavy species in the upper atmosphere of HD 209458b 


Instrument Species 

Species 

Species 

Species Ref. 

mode feature (A) 

feature (A) 

feature (A) 

feature (A) 


Absorption depths from transit light curve 

FUV observations 



OI 1304 

CII 1335 

Si III 1206.5 

Si IV 1393.7 



triplet 

doublet 

line 

line 


STIS G140L 

12.8±4.5% 

7.5±3.5% 

0.0/qq% 

0.0/®;®% 

Vidal-Madjar et al. (2004)“ 

STIS G140L 

10.5±4.4% 

7.4±4.7% 

confirmed 

confirmed 

Ben-Jaffel &; Hosseini (2010 ) a 

STIS G140M 

n/a 

n/a 

1.7±18.7% 

n/a 

this work a 


Invalid depth estimates from COS 2009 dataset from in-transit data against data from other dates 



CII 1335 

Si III 1206.5 

Si IV 1393.7 



doublet 

line 

line 


COS G130M 

7.8±1.3% 

8.2±1.4% 

0.00±0.01% 

Linsky et al. (2010 ) b 

", this work 

9.3±2.4% 

8.6±2.7% 

3.1±4.2% 

vs off-tr. ave, photon noise c 

", this work 

5.3±2.9% 

-1.9±3.5% 

1.5±4.9% 

vs <f> = 0.27, photon noise 6 * 

", this work 

9.3±8.4% 

8.6±9.5% 

3.1±12.1% 

vs off-tr. ave, w intra-visit var. e 

", this work 

5.3±9.8% 

-1.9±12.0% 

1.5±16.8% 

vs (f) = 0.27, w intra-visit var. e 


“ The STIS G140L and G140M results were respectively derived from 4 and 3 transit observations. 

b Invalid ’’transit depth” reported bv ILinskv et a DOT), not from a transit light curve but from the average of data with the planet 
off-transit at orbital phases <j> = 0.27, 0.72 and 0.49, and with errors for photon noise only. 

c As in item 6, for comparison with the Linsky et al. results, with photon noise only. 

d Sample invalid ”transit depth” from contrasting the in-transit data against off-transit observation at cf) = 0.27. A negative detection 
is obtained in Si III, and a smaller transit depth is obtained in CII. 

e Invalid ” transit depths” derived as in items c or d but including the intra-visit standard-deviation variations (Table 2) in the errors. 
No previous HD 209458b UV transit study has included this variation. Full transit light curves are needed for a proper assessment 
of these variations and thus of the true potential of the COS transit measurements of HD 209458b. 


the upper atmosphere derive from dissociation products 
in the lower atmosphere, issues of metallicity and C/O 
ratio are directly testable with FUV observations - pro¬ 
vided high quality data resolving the transit absorption 
line shapes are obtained. 

The presence of Si III in the upper atmosphere of 
HD 209458b would have implied that silicon gas is not ef¬ 
fectively d epleted by cloud conden s ation in the lower at - 
mosphere (IKoskinen et al.ll2013albt ILavvas et al.1 l2014h . 
This turns out to possibly still be the case as demon¬ 
strated by the recent detection of 6.2 ± 2.9% absorp¬ 
tion (or 8.8 ± 2.1% if there was still post-egress ab¬ 
sorption) by extended Mg I on this_planet at near-solar 
abundance (|Vidal-Madiar et aLll2003lh because magne¬ 
sium and silicon condense together int o si li cate g rains 
of fosterite and enstatite (lYisscher. Lodders fo Feglevl 
l201d h The detection of magnesium in the thermo¬ 
sphere of HD 209458b indicates that the balance of day- 
to-ni g ht temperature (Showman et al . 2009; Mose s et al.l 
120111: iCrossfield et al.ll2012fl aga inst the strong verti¬ 
cal and horizontal 3D dynamics dShowman et al.l [20091 : 
iSoiegel et al.ll2009t iParmentier et al.l 120131) impedes sig¬ 
nificant condensation and settling of refractory si li¬ 
cate species on t his planet dKoskinen et al.1 l2013allbl : 
ILavvas et al.l [2014b If condensation takes place, such 
as on the nightside, the transport of fine silicate grains 
onto hotter high-altitude dayside regions may allow for 
effective sublimation on this hot Jupiter. Such a strong 
transport of fine grains has been found to be possible 
on HD 209458b based on 3D global circulation modeling 
that included test particles and found a strong circula¬ 
tion with large effective eddy coefficients of fo r exam ple 
K zz ~ 10 lo cm 2 s _1 at P ~ 1 mbar dParmentier et all 
120131) . 

With respect to the FUV and NUV observations and 
the upper-atmosphere modeling, the disagreement lies in 


the predicted Si III and observed Mg I ionization states, 
such that silicon might be present mainly as either Si I 
or Si II. In the latter case charge-exchange reaction rates 
invoked to convert part of the Si II into Si III would need 
to have been overestimated so that no significant Si III 
1206.5 A absorption is detected in transit. The charge ex¬ 
change of Si II with protons ( Si + + p + —> Si ++ + H) was 
first proposed bv ILinskv et al.l (120101) as the main source 
of Si III since this reaction has been found to be relevant 
for the Sun. In their detailed upper-atmosphere mod- 
els for HD 209458b th at included heavy metal species, 
iKoskinen et al.l (l2013al lbl) confirmed that this reaction is 
key in producing detectable amounts of Si III to fit the 
~8% transit depths reported by Linsky et al. from the 
COS dataset. However, our new STIS G140M negative 
finding for Si III on this planet that is consistent w ith 
the_previous negative detection by iVidal-Madiar et ahl 
(l200l (Table 5), poses new questions into the model¬ 
ing and/or the state of the upper atmosphere of this 
planet. It may be that the model temperatures are too 
high since temperatures above ^15,000 K are required 
for the Sun to produce significant Si III abundances 
(iBaliunas fc Butlerill980 ) and these temperatures do not 
necessarily apply to this planet. The model tempera¬ 
ture depends on the chemical composition, on the stellar 
XUV input and the heating efficiency, and on the velocity 
of the outflow since it can adiabatically cool the atmo¬ 
sphere, although there are always uncertainties in these 
calculations. We find in the literature that only models 
with high stellar XUV input and high heating efficiency 
seem to reach high temperat ures o f 15,000 K and beyond, 
such as all the model cases inlTian et al] (|2005f) , the sub¬ 
sonic cases in iGarcfa Muhoi ( 200711 . and the models in 
IKoskinen et al.l (l2013aU bli~ that explicitly assume a xlO 
or xlOO higher solar XUV input for the star. Other 
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models reac h about 12,0 00-14,000 K peak temperatures, 
such as the lYellel (12004T) mode l, the super-sonic cases in 
IGarcia Munoz] (|2007D . and the IKoskinen et aD (l2013al fbl) 
models with solar XUV input. In the modeling of the 
upper atmosphere of a non-magnetized HD 2094 58b and 
its interaction with the stellar wind, iMurrav-Clav et ahl 
(120091) find that the peak temperature reaches ~10,000 K 
for a solar XUV flux, while for a 1000 times larger XUV 
flux typical of young T-Tauri stars the peak tempera¬ 
ture reaches only slightly larger values due to radiative 
HI Lyman-a cooling rather than by the expansion of the 
gas. In a separate work we are finding evidence that 
the star HD 209458 has a surprising long-term low ac¬ 
tivity, and this indicates that the star provides a sig¬ 
nificantly lower XUV inp ut to t he planet than assumed 
in all previous modeling (jBallest.er et, al.ll2015D . In that 
work we apply the detailed photochemical h ot-Ju piter 
upper-atmosphere models of IGarcia Munozs (|2007l) that 
include detailed photochemistry for both ion and neu¬ 
tral species of H, D, He, O, C and N and for which we 
have updated our best current estimate of the stellar X- 
ray and UV input. The modeling finds indeed peak at¬ 
mospheric temperatures (for solar abundances) of only 
^10,000 K. Given the Mg I detect ion, a nd acc ompanying 
non-detection of Mg II (lVidal-Madiar et all 120131) , the 
dominant silicon ionization state may be Si I since Mg I 
and Si I have similar ionization potentials. Nevertheless, 
the finding that Mg I dominates in the upper atmosphere 
of this planet is highly unexpected, since significant pho¬ 
toionization (A < 1621 A) should be present while Mg II 
was not positively detected. This seems to require very 
large electron densities for the dielectric recombination 
required to counteract photo-ionization: f or ex ample, 
at the reference altitude of r — 3R P iBourrier et al.l 
(120141) require at least two orders of magnitude higher- 
electron densities than estimated from current hydro¬ 
dynam ic models (IGarcia M unod [20071: iKoskinen^.al. 
2013allbl : iGuol 120131: lLavvas et al.l 12014 iBallester et al. 
20151 ). This is something that is rather hard to explain. 

Another finding is that the Mg I 2853.0 A line ab¬ 
sorption seems to b e blue-shifted to -62 to -19kms _1 
(lVidal-Madiar et aLll2013D . This finding is quite com¬ 
pelling to some and unexpected to others, as per the 
discussion in the Introduction related to the HI Lyman- 
a line. Again, the large electron densities required for 
dielectronic recombination and the fast outflows needed 
at the exobase to explain the Mg I bl ue-s hif ted a bsorp- 
tion by stellar radiation pressure (IBourrier et al.l 120141) 
are difficult to reconcile with current upper atmosphere 
models. Yet, if confirmed, the Mg I blue-slrifted absorp¬ 
tion reveals a net anti-solar motion that may apply to one 
or more components of the upper atmosphere. Do the OI 
and CII species also show a net blue-shifted absorption? 
Are t he motions of the ions and neutrals fully coupled 
(e.g., IGarcia Mufiozl 120071 iTrammell et al.l[2011b lAdamsI 

120111) or are they decouple d at large enough dis tances 
and low collision rates (e.g.. IKoskinen et al.ll2014l )? The 
latter could be expected from effects by different radi¬ 
ation pressure, charge exchange with the stellar wind 
protons, close and open planetary magnetic-field lines 
and decaying field strength with distance, and net mag- 
netospheric currents and convection. Velocity-resolved 
transit data on the CII 1335 and OI 1304 A lines, even 


if only separating the blue- and red-component absorp¬ 
tions, would provide much needed tools for characterizing 
the state of the upper atmosphere of HD 209458b that is 
still an enigma. Such data can be obtained with COS. 

5. CONCLUSIONS 

We have presented new STIS G140M transit data re¬ 
inforces the original results with STIS G140L that Si III 
ions are not detected in the thermosphere of HD 209458b. 
Silicon may still be present in the upper atmosphere, but 
at a lower ionization state that may even be Si I, based 
on the independent HST detection of extended Mg I on 
this planet which also indicates that silicate condensation 
and cold trapping is not 100% effective in t he lower atmo¬ 
sphere. The models for HD 209458b by Koskinc n et al.1 
(|2013allbh that include refractory species already predict 
Si II to be dominant silicon species in the upper atmo¬ 
sphere, so what is needed as a first step is to limit the 
yield for Si III. A second step much harder to address 
would be whether the silicon is in the Si II or Si I state, 
given the positi ve Mg I but negativ e Mg II detection on 
HD 209458b by lVidal-Madiar et ahl (120131 ) and that Si I 
and Mg I have similar ionization potentials. The lack of 
a positive Si III detection in the extended atmosphere of 
the planet reported in this work, together with the pos¬ 
itive Mg I detection and the long -ter m low activity now 
being identified for the star IBallester et al.l (120151) indi¬ 
cate that a revision to our present understanding of the 
upper atmosphere of this planet is needed. Much new 
insight can be obtained with new COS FUV transit ob¬ 
servations. 

The existing COS G130M observations of the 
HD 209458 system were not standard transit observa¬ 
tions and significant stellar flux variations found in 
the data invalidate previously reported transit results. 
Proper FUV transit observations with COS are needed, 
sampling the transit light curve (at a single grating wave¬ 
length setting). These observations will independently 
confirm and measure transit depths in the OI and CII 
species at high S/N, and furthermore resolve the velocity 
distribution of the OI and CII species with major impli¬ 
cations on upper atmospheric dynamics, energetics and 
magnetosplieric processes, as well as on the species abun¬ 
dances that are also relevant for the characterization of 
the planet and its lower atmosphere. In this frame, the 
detection of OI and C II in the upper atmosphere of 
HD 209458b using only unresolved lines already requires 
solar or supersolar abundances that are in disagreement 
with a very low metallicity as recently considered in the 
literature. In the near future, upper atmospheric UV 
studies of exoplanets may prove key for unraveling the 
so far elusive properties of low mass exoplanets. 
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